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We analytically study heat conduction in a chain with interparticle interaction V{x) = A[l — 
cos(j:)] and harmonic on-site potential. We start with each site of the system connected to a 
Langevin heat bath, and investigate the case of small coupling for the interior sites in order to 
understand the behavior of the system with thermal reservoirs at the boundaries only. We study, 
in a perturbative analysis, the heat current in the steady state of the one-dimensional system with 
weak interparticle potential. We obtain an expression for the thermal conductivity, compare the low 
and high temperature regimes, and show that, as we turn off the couplings with the interior heat 
baths, there is a "phase transition:" the Fourier's law holds only at high temperatures. 
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The understanding of heat conduction in a lattice sys- 
tem of interacting particles has become a challenging 

Hoblem of statistical physics, even in the 11? context 
. A central issue is finding a model Hamiltonian sys- 
tem for which Fourier's law holds. One of the first works 
in this subject was the (rigorous) study of the harmonic 
chain of interacting oscillators coupled to heat baths at 
the boundaries The authors show that the heat cur- 
rent is independent of the length of the chain, i.e., the 
Fourier's law does not hold. Since then, many (very often 
conflicting) works have been devoted to the problem, in 
particular, to investigations on the effects of nonlinearity 
and external potentials in the behavior of the heat cur- 
rent. We recall some results. In the authors show 
that the conductivity is anomalous (i.e., it diverges) in 
any one-dimensional momentum conserving system, but 
in and a momentum conserving system with fi- 
nite conductivity is presented. In 0, the authors claim 
that the anharmonicity of the on-site potential is a suf- 
ficient condition for a finite thermal conductivity, but 
in 0, it is shown to be wrong. Almost all the results 
are obtained by means of computer simulations, and, as 
emphasized in j^, besides the difficulty to arrive at cor- 
rect conclusions from numerical studies, several works 
use the Green-Kubo formula for the conductivity, a for- 
mula which has never been rigorously established for this 
context. In short, more accurate studies are necessary. 

In this scenario, the harmonic Hamiltonian chain of 
oscillators has been revisited quite recently but for 
the case of each site connected to a thermal reservoir. 
The steady state is rigorously computed in the "self- 
consistent" condition, which means with no heat flow 
between an inner site and its reservoir. In such a model, 
the Fourier's law holds. 

The present paper is addressed to the following is- 
sues: (i) the development of new analytical methods of 
modeling the heat conduction problem; (ii) the search 
for a system with normal conductivity and with, say, a 
small anharmonic potential (the problem which inspired 
the first investigation of Fermi, Pasta and Ulam [HI); 



(iii) the understanding of the temperature role on the 
thermal conductivity of chains with soft anharmonicity 
such as V = 1 — cos((;i+i — Qj) (there is a recent de- 
bate, with positions against jlj and in favor of a 
phase transition in the rotor model - i.e. finite thermal 
conductivity for large T, and infinite one for small T). 
Here, we extend the approach and techniques previously 
developed in in order to treat a chain with ther- 
mal reservoirs at the boundaries only: now we consider 
different coupling constants among reservoirs and sites, 
and investigate the limit of the coupling with the interior 
heat bath taken to zero. Our approach is quite general, 
but we focus on the case of a chain of oscillators with a 
harmonic on-site potential and interparticle interaction 
V = A[l — cos(gi+i — Qi)]. We obtain (in a perturba- 
tive analysis) an expression for the thermal conductivity 
and investigate the Fourier's law: for our model, as we 
turn off the the couplings between inner sites and their 
reservoirs, it holds only at high temperatures. 

Now we introduce the model. We consider the 
Langevin dynamics of an anharmonic crystal with 
stochastic heat bath at each site. Precisely, we start from 
N oscillators with Hamiltonian 

N N 

(1) 

(d = 1 and next-neighbor interactions are assumed later) 
where M > 0, with time evolution, for j = 1, . . . , TV , 

OH 1/2 

dqj = pjdt; dpj = -- — dt - CjPjdt + dBf, (2) 

OQj 

where Bj are independent Wiener processes; Q is the 
heat bath coupling for the j*'* site; and jj = 2QTj, where 
Tj is the temperature of the j*'' heat bath. 

As usual, we define the energy of the oscillator j as 

Hj(.q,p) = \p] + U^'Hqj) + ^^'^fe - qi), (3) 
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where the expression for U^^'' and U^'^^ follow immedi- 
ately from (Q) and X^^i — ^- Then, we get 

where (•) denotes the expectation with respect to the 
noise distribution, and 

{RM=Q{T,-{p'j)) (5) 

gives the energy flux from the j*'' reservoir to the j*^ site. 
The remaining terms are related to the energy current 
inside the system and they are given by 

^,^ = ;^vc/(^)fe-90^; (6) 
i>j ^ 

describes the heat flow from the j*'' to the Z"* sites; 

is obtained from the formula for by changing 
I with 2- It is useful to introduce the phase-space vector 
(\) — {q,p) with 2N coordinates and write the equation 
for the dynamics (2) as 

^ = -A0-[/(2)' + cr77, (7) 
where A and a are 2N x 2N matrices given by 

/ above is the unit N x N matrix; and A4,T,T are di- 
agonal N X N matrices: Aiji = MSji , Tji — Cj^ji j 
Tji — TjSji. r] are independent white-noises; f/'^'' is the 
derivative of the J/^^^ term in H in relation to q (note 
that its contribution to (pk is nonzero only for k > TV). 

To study the dynamics we adopt the following strat- 
egy. First, we consider the system with J/^^-* — 0, and 
stay with N independent sites connected, each one, to 
a heat bath. To recover the original dynamical system, 
we introduce the interaction among the sites and calcu- 
late the changes using techniques of stochastic differential 
equations. The solution of O above with J/'^' = 0, is 
the Ornstein-Uhlenbeck process 

Jo 

For simplicity we take (/)(0) = 0. The covariance of this 
Gaussian process is 

C{t,t) = / ds e-'^a^e-'^^. 
Jo 



From an easy computation (e.g. diagonalizing A), it fol- 
lows that (for a single site (pj) 

exp {-tA) = e^*^ cosh{tpj) | ^ g ^ ^ 

-^(14))' <'°' 

1/2 

Pj = ((Cj/2)2 — M) ; the expressions for involv- 
ing 2N X 2N matrices are immediate. We assume that 
Cj/2, M > 0. If (0/2)^ > Af, then is real; 
otherwise, pj is pure imaginary, but it does not spoil 
the dynamics: cosh(tp) in the formula above becomes 
cosh(t • ip') — cos(tp'), etc. In this case (i.e., with 
= 0), as t ^ oo we have a convergence to equi- 
librium and the stationary state is Gaussian, with mean 
zero and covariance 

C = l°^ds e-^a^e-^" ^ I ^ ), (11) 

where T is a diagonal matrix with elements TiSij. To 
introduce the anharmonic coupling potential, we use the 
Girsanov theorem |l4| , which establishes a measure v for 
the complete process as an integral representation in- 
volving the measure pc associated to the process without 
the potential C/^^-* . Precisely, for any measurable set R, it 
states that p{R) — EQ{lfjZ{t)), where Eq is the expecta- 
tion for PC (the measure for the process with 1/*^^^ — 0); 
Ir denotes the characteristic function, and 

Z{t) — exp ^/J u • — i /g u'^ds^ ; 

7y'u, = -V,_A,C/(2); (12) 

the inner products above are in and Vfc means the 
derivative in relation to (pk- From ||HJ| and the expression 
above for Ui, we have that Ui is nonvanishing only for i > 
N (i.e., i G [N+1, 7V-f 2, . . . , 2N]). In what follows we will 
use the index notation: i for index values in the set [N + 
1, iV-l-2, . . . , 2iV], j for values in the set [1, 2, . . . , TV], and 
k for values in [1,2,..., 2iV]. We will also be restricted 
to next-neighbor interactions, i.e., we take 

[/(2) ^ l;7(2)(0i - ^2) + ^c/(2)(<^^_i - 0^) + 

N-l 

where U'^'^\x) = A[l — cos(x)]. Using the A: index 
notation above, we may rewrite the stochastic equations 
for the decoupled process (where C/'^^ = 0) as 

1 /2 

d(f)j = -Ajkipkdt; d(f)i = -Aik(j)kdt + -f^ dB^; (13) 

the sum over k (in [1, 2 . . . , 2iV]) is assumed above (as 
well as obvious sum over some indices in what follows). 



Now, let us make explicit the terms in Z{t). We have 

where we dropped out the upper index in t/^^^ above 
(and in what follows); U-_j^ means the derivative in re- 
lation to (j)i-N, i-e., U'j = U'{(j)j — (jij+i) — U'{(hj-i_ — 
U'{(f>j) — — Asin(0j). From Ito formula it 
follows that 



of the second term in the r.h.s of the equation above dis- 
appear in the correlation formula in the limit of i — > oo 
(recall that we must take this limit i ^ oo in order to 



reach the steady state). Writing Z{t) 
vious formula becomes 



(U' ~U' ) ^ 



jdt (jjll 



where C/' 



U'{4ij — etc; and j — i ~ N in the 

expressions above. Hence, we get 

~ ;7j_i JA,fe0fe(s)ds) , (14) 



and, for the term. 



ds = -\- \ [C/j^^.+i - C/j_iJ^(s)ds, (15) 



where, again, j — i — N an the sum over i (and 
so j) is assumed above. Then, for the correlation 
functions, we obtain an integral representation involv- 
ing a "perturbative" potential and a Gaussian measure. 
E.g., for the two-point function we get {(f>k(ti)4>q{t2)) = 
Af J Mti)Mt2)Z{t)dnci(j)), tiM < t- where Z{t) = 
e~^, W is described by the several terms presented by 
the expressions above; A/" is the normalization. 

The heat flow in the steady state is related to the for- 
mula jnjl. Precisely, for the case of next-neighbor inter- 
actions, the average over the stationary distribution for 
the current T^^v+i is obtained as the limit 

limt_oo(-?\.^t,+i) = limt^oo 
![K,,+,{ct>u + c^u+i)/2mZ{t)dtic{<ty)/ j Z{t)dt,c{^) , 

where v = u — N {u > N, obviously). Now, we 
carry out the computation. First note that C{t, s), given 
by (fTn|l - l(TT|l . may be written as (for t > s) C{t,s) — 
exp(-(t - s)A)C + O (exp[-(t -I- s)C/2]), and the effects 



-\w 



, the pre- 



lim {Ty^y + l) 
t—>-oo 



lim 

t — >OQ 



-xw 



d^icicl)) 



il given by the product of U' and (j) described above. Up 
to first order in A (i.e., for weak interaction between two 
sites), we have {J^v^v+i) = {^)c + —^W)c, where 
' {■)c means the average in respect to d^c; (■;')c means 
the truncated expectation. It is easy to see that (ri)c 
vanishes: U' depends on or (pv+i, and {4>v{t)4>u{t))c = 
Cy^u{t,t) = (for any v < N and u > N). The terms 
in ~XW are given by those describing UidBi ifTHl . dis- 
carding that one involving (/)(0) which vanishes in the 
computation as t — > oo (C(i, 0) — > as i ^ oo). Note 
that we stay involved with expressions such as 

(C^;.+i0u(t);C/i^,+i0j(s))c = 

_i ^gi0„+i(t)e-»<A-(*)0„(t);e*'>^+i(")e-''^^(''Vj(s)) + • . ■ 
i.e., with integrals hke / e*(''i+ --+'"''''^c?^c ~ 

Je'{hi + h2+h,)-4>^^^J^r{h,+h,+hey^^^^^ (0„(i) is 

obtained from e*'' "^, obviously, by taking the derivative 
in relation to hu(t) and making h = 0; the same for 
(f)~j{s)). After these integrations in cj), we get expressions 
such as 



lim / ds C„j+i(t, s)C , , :;(t, s)exp[- 



a,l3 



Ca,/3(iQ, ^/3)/2], 



where ta,t/3 are t or s. For the case of small tempera- 
tures Tj, we have exppc^/j] « 1. Thus, after all (j) and s 
integrations (taking also the limit t —^ oo), we obtain 



-A2 



2M(C„ + C.+i) 



Cv rri Cv + 1 rji 
J'U + l T Iv 



+ 1 



Cv 



(16) 



(to simplify the notation we write Ty—^y+i in- 
stead of \mit~>oo{J^v~>v+i))- For the high tempera- 
ture regime, we have to deal with expressions like 
limt_oo/o /(i,s)e-^9(*'")£is, where 6* cx T. We use the 
Laplace method 0| to get their asymptotic behavior 
9 Qo. We obtain 



4 



with slight changes for the terms with v = 1 and v — N. 

The steady state is characterized by (dHj/dt) = 
0. Using this expression and also that hnit^oo('^f) = 
Ti-N (taking the dominant contribution), which gives 
limt^oo {Rj{t)) = for the interior sites j (i.e. the "self- 
consistent condition"), we have T\^i = = ... = 
Tn-i^n = ^- Hence, for Qj+i — C,j small, summing up 
!Fi^2 + -^2^3 + ... we obtain, for small temperatures, 

T ■ 2M(Ci + 2C2 + • • • + Kn-i + Cn) ~ A2(Ti - Tjv). 

For uniform C, we have the Fourier's law 

T^x{Ti-Tn)/{N-1), X-AV4CM. 

As we make the inner couplings smaller and smaller we 
lose the factor N (which comes from C2 + • • ■ + Cjv-i) and 
the Fourier's law does not hold anymore. For the case 
of high temperatures the sum of all JF^^t,+i gives us, for 
Cv+i - (y small. 



A2 



-Tn/M 



N 



-Ti/M 



which (essentially) does not depend on inner heat bath 
couplings, and so, the Fourier's law still holds when we 
make them smaller and smaller. Taking (^n — Ci — Cy 
and Tn — Ti + d {d small) the expression above becomes 



4MC 



'T/M 



(Ti - Tn) 
{N-1) ■ 



(18) 



where e'^/^^ = -[e-iT.+S)/M _ e-T'i/Af]/[(Ti + 5)/M - 
Ti/M], i.e., the conductivity decays exponentially at high 
temperatures, as in the rotor model Thus, still con- 
cerning the chain of rotators and the recent debate about 
the existence of a phase transition , j based on our 
results with cosine interactions, we believe in a divergent 
conductivity in low temperatures, as claimed in |l2l |. 

To argue about the reliability of our treatment, we 
recall some previous related works where the perturba- 
tive analysis gives the same result of the rigorous treat- 
ment. For the simpler case of the harmonic chain of 
oscillators with a bath at each site and identical next- 
neighbor interactions, a first order perturbative analy- 
sis (for weak interactions, in a similar approach to 
that described here) gives the same result of the com- 
plete and rigorous treatment And following the pro- 
cedures described here (analyzing different Qj), one may 
see that the perturbative result for this harmonic chain, 
in the limit of zero coupling between reservoirs and in- 
ner sites, will lead to the rigorous result obtained for 
the harmonic chain with thermal baths at the bound- 
aries Q . We still recall some previous works considering 
nonconservative stochastic Langevin systems (in contact 
with thermal reservoirs at the same temperature), but in- 
volving similar integral expressions for the correlations. 



There, the time decay of the two and/or four-point func- 
tions is detailed investigated in the regions of low and 
high temperatures. For the low temperature regime and 
weak interaction among the sites, we rigorously prove jl^ 
that the complete treatment of the two and four-point 
functions adds only small corrections to the perturbative 
results 0] . For the same nonconservative system at high 
temperature, we developed a cluster expansion [l^ which 
supports the perturbative analysis [l9j . 

In short, in the present letter we develop an analytical 
method of modeling the heat conduction problem and 
study the heat current at the steady state of an anhar- 
monic chain with weak interparticle (cosine) potential in 
order to investigate an old problem of heat conduction: 
may small anharmonic interactions lead to normal con- 
ductivity? We show that, if we keep the thermal reser- 
voirs at the boundaries only, the Fourier's law holds for 
high but not for low temperatures. 
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